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The charge transport of a serially coupled quantum dots (SCQD) connected to the metallic 
electrodes is theoretically investigated in the Coulomb blockade regime. A closed-form expression 
for the tunneling current of SCQD in the weak interdot hopping limit is obtained by solving an 
extended two-site Hubbard model via the Green's function method. We use this expression to 
investigate spin current rectification, negative differential conductance, and coherent tunneling in 
the nonlinear response regime. The current rectification arising from the space symmetry breaking 
of SCQD is suppressed by increasing temperature. The calculation of SCQD is extended to the case 
of muhiple parallel SCQDs for studying the charge ratchet effect and SCQD with muhiple levels. 
In the hnear response regime, the functionaUties of spin filter and low-temperature current filter are 
demonstrated to coexist in this system. It is further demonstrated that two-electron spin singlet 
and triplet states can be readily resolved from the measurement of Seebeck coefficient rather than 
that of electrical conductance. 



I. INTRODUCTION 

Tunneling current through individual quantum dots 
(QDs) with discrete levels exhibits a kaleidoscope of in- 
teresting physics such as the Kondo effect, Fano reso- 
nance, and Coulomb blockade.^ Recently, a serially cou- 
pled quantum-dot (SCQD) system (the simplest artifi- 
cial molecule) was proposed as a spin filter based on 
the spin-blockade process for appUcation in spintron- 
ics and quantum computing. The transport proper- 
ties of the SCQD in the Coulomb blockade regime in- 
clude current rectification,^"^ negative differential con- 
ductance (NDC),^~^ nonthermal broadening of electri- 
cal conductance,^'^ and coherent tunneling for the SCQD 
with degenerate energy levels (quantum dot "helium").^ 
Many theoretical efforts have been devoted to studying 
them.'^"-'^* Nevertheless, there still lacks a comprehen- 
sive theory to explain these phenomena in a systematic 
way.-^^ The early studies of SCQD focused on the co- 
herent transport behavior without spin dependence. ^'^'^ 
For the application of spintronics and experimental ob- 
servations of spin-dependent tunneling current, the cur- 
rent rectification arising from the Pauli spin blockade and 
NDC were theoretically studied by several groups. ^^"^"^ 
The transport properties of SCQD embedded in a ma- 
trix connected to ferromagnetic electrodes were also 
studied. However, theoretical study of the nonther- 
mal broadening of electrical conductance of SCQD in the 
presence of electron Coulomb interactions has not been 
reported. 

Using the Green's function (GF) technique we have 
solved an extended Hubbard model, which includes 
the interdot Coulomb interactions as well as intradot 
Coulomb interactions for a coupled quantum dot system, 
in the weak interdot hopping strength limit. The de- 
rived closed-form solution for the transmission factor has 
eight spin-charge configurations and 16 resonant chan- 



nels. With this theory we can provide quantitative anal- 
ysis for current rectification arising from coherent tun- 
neling with spin blockade, NDC, and nonthermal broad- 
ening effect of tunneling current resulting from the off- 
resonant energy levels. All our predictions are in good 
agreement with available experimental data.^~^ In ad- 
dition, we demonstrate that the SCQD junction system 
can be used as a low-temperature current filter and spin 
filter simultaneously. 

To depict the charge ratchet effect arising from other 
QDs surrounding the SCQD^° and the case of SCQD 
with multiple levels (such as the Si-based SCQD, where 
the multi-valleyed nature of Si gives rise to closely-spaced 
energy levels in each Si QD,^^'^^, we extend our calcu- 
lation to a four-level Anderson model by using a sim- 
ple and physical picture. Unlike GaAs QD system with 
large nuclear and electron spin interactions,^ Si QDs have 
longer electron spin coherent time, which is an impor- 
tant consideration for cjuantum computing. Owing to 
traps and impurities of Si QDs and their multi-valleyed 
characteristics^^'^^, we need to consider a Hamiltonian 
beyond the two-level Anderson model. However, based 
on the conventional theoretical framework,^ it is compli- 
cated to deal with carrier transport properties of multiple 
QDs. In this paper, we propose a new and comparatively 
easy approach, which allows us to derive a closed form 
solution of tunneling current for arbitrary QD numbers 
in the limit of weak interdot hopping strength. Based on 
this closed-form expression of tunneling current, we can 
fully analyze the tunneling current spectrum of multiple- 
SCQD system. 

Recently, many theoretical efforts have been devoted to 
the studies of thermoelectric properties of QDs in quest of 
highly efficient thermoelectrical materials. ^'^^^^ However, 
these studies have focused only on the thermoelectrical 
properties of a single QD^'^'^^ or parallel QDs.^^'^'' To re- 
duce the temperature gradient across the QD junction. 
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it is essential to study the case of N coupled QDs in se- 
rial. In this paper, we also investigate the thermoelectric 
effect of SCQD. We find that the Seebeck coefficient is 
much larger in the spatially symmetric SCQD than the 
spatially asymmetric SCQD. Furthermore, we find that 
the measurement of Seebeck coefficient provides an alter- 
native means to distinguish the spin singlet and triplet 
states in the linear response regime. 

II. FORMALISM 

The inset of Fig. I illustrates the model system of con- 
cern, showing a serially coupled quantum dots connected 
to metallic electrodes. The electron Hamiltonian can be 
described by a two-level Anderson model:^ 

k,cT k.cT £,a 

+ ^Utd\^^di,^d\^^di^^ + ]^ ^ 

+ ^ \V^.Aal. ,^dA,a + Vi.Brfs.^&k.a] + h.C. 
k,<T 

+ ^tl2(4,.cis,a-f (1) 

where the first two terms describe, respectively, the 
free electron gas of the left and right metallic electrodes. 
The operators in the system Hamiltonian are defined as: 
'^k cr ('^k,cr) creates (destroys) an electron of momentum 
k and spin a with energy et in the left metallic elec- 
trode. ^ (^k.o-) creates (destroys) an electron in the 
right electrode. d\ ^ (d^.a) creates (destroys) an electron 
in the £th. dot (dot A or dot B). Ui and Uij describe the 
Coulomb interactions inside the £th dot and between the 
^th and jth dots, respectively. For small semiconductor 
QDs (with size ~nm) , the orbital energy level separation 
of individual QD is much larger than and thermal en- 
ergy ksT. This allows us to consider only one energy 
level for each dot in Eq. (1). describes the cou- 

pling between the band states of electrodes and state of 
dot A (B) with energy level Ei {E^)- The last two terms 
describe the electron hopping between two dots. 

Using the Keldysh-Green's function technique, ^'^ we 
obtain the tunneling current expression of SCQD (see 
Appendix) as 

J = ^ 1 &r(6)[/L(6) -/«(£)], (2) 

where T(e) = T L{e)V R{e){Ai2 + ^2i)/2 is the transmis- 
sion factor. T(^L^}i{e) denote the tunnel rates from the 
left electrode to dot A and from the right electrode to dot 
B, respectively. fL{R){f) = l/[e("-''^(«))/''«'^^(«' + 1] de- 
notes the Fermi distribution function for the left (right) 
electrode. The chemical potential difference between 
these two electrodes is related to ji^ — I^r = eAVa. 
denotes the equilibrium temperature of the left (right) 



electrode, e and h denote the electron charge and Plank's 
constant, respectively. For simplicity, we consider the 
wide- band limit: r^(e) = Tg. The "resonant function" 
v4f J- of the transmission factor can be calculated by the 
on-site retarded Green's function {G\ ^{e)) and the lesser 
Green's function {Gf^{e)) (see Appendix). 

After straightforward algebra, we obtain an expression 
for the resonant function 
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Aj(e) = t?2 E P^n/\^.n\^; {i ^ j), (3) 
ni—1 

where the denominators for the eight configurations are: 
(i) Hi = Hi^J-j — with both dots empty, (ii) 112 = 
{^i — Uij) (/ij — Uj)— , with dot £ empty and dot j filled 
by one electron with spin ct, (iii) Ha = (/i^ — t/fj)(/ij — 
Ujj)—ti2 with dot £ empty and dot j filled by one electron 
with spin cr, (iv) 114 = {^J^t. — '^Uij)(^j — Uj — Uj^i) — 
with dot ^ is empty and dot j filled by two electrons, 
(v) 115 = (/^£ - Ut){nj - Ujj) - with dot j empty 
and dot £ filled by one electron with spin a, (vi) Ifg = 
{fii — Ui — U£j){fij — Uj — Ujj) — with both dots 
filled by one electron with spin ct, (vii) Ily = {^i — Ui — 
Uij){^j — 2Uj,i) — with dot i filled by one electron 
with spin a and dot j filled by one electron with spin cr, 
and (viii) ITg = {^i — Ue — 2Uij){^j,j—Uj — 2Uj,i)~tl2 with 
dot i filled by one electron with spin a and dot j filled 
by two electrons, fie = e — Ei + iTe/2. The numerators 
PrnS denote the probability factors for various sin-charge 
configurations. They are pi = (1 — iV^_g.)(l — iVj^cr^ A^j,ff + 
Cj),p2 = (l-iV£,ff)(iV,,ff-c,),p3 = (i-iVf,ff)()Vj,a-Cj), 

P4 = (1 - Ni^ff)Cj, P5 Ni^^{l - Nj_„ - Nj_g + Cj), P6 = 

Ne,&{Nj^g- - Cj), pj = Ni^g{Nj,„ ~ Cj), and ps = A^^.^c^ 
((T denotes the opposite of cr), where Ng^^ and ct denote 
the thermally averaged one-particle occupation number 
and two-particle correlation function, respectively. They 
can be obtained by 

N,,. = M = I ^G</e), (4) 

Q = / ^Gffie), (5) 

where Gfg{t) and Gl'f{e) denote, respectively, the on- 
site one-particle and two-particle lesser Green's func- 
tions, which can be calculated by the equation-of-motion 
method. 

The expression of Eq. (3) is valid in the Coulomb 
blockade regime, but not in the Kondo regime, be- 
cause we did not take into account the effect of electron 
Coulomb interaction on the tunneling rate (F^) which 
arises from the coupling between the electrodes and the 
QD. Eq. (3), also valid in the limit of ti2/U ^ 1, is cor- 
rect up to the second order in ti2. This is sufficient for 
our analysis of weakly coupled dot, since in the SCQD 
system the value of hopping strength ti2 is much smaller 
than all the other energy scales. When ii2/r£ ^ I, we 
can safely neglect in the lesser Green's functions the cor- 
rections coming from the neighboring dot. Eqs. (4) and 
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(5) can then be rewritten in terms of the on-site retarded 
Green's functions: Ng^^ = — (I/tt) / c?e/f (e)ImG^ ^(e), 

and C£ = —(I/tt) J de/£(e)ImG^'^(e), where the retarded 
Green's functions G\ g{e) and G^'^(e) are given in Eqs. 
(A. 16) and (A. 17). They arc solved self-consistently. We 
find that the resonant channels of Ai^j arc related to the 
off-diagonal one-particle Green's function, which is given 
by 



Pi 



(6) 



P2 



+ ?1 

+ El 

{fii - 2Ue.j){^ij - Uj - Uj^e) - i?2 
+ El 

+ El 

+ El 

+ El 

(/Xf -Ut- 2C/£j)(At, - Uj - 2Uj,i) - ■ 

To study thermoelectric properties in the linear re- 
sponse regime, we rewrite Eq. (2) as 



J = CnAV + C12AT, 



(7) 



where AT — Tl — Tn > is the temperature difference 
between two electrodes. Coefficients in Eq. (7) are given 
by 



2e2 f 
^11 = — y der{e){ 

2e f .(9f{e) 
C12 = deT{e){-^)E, 



dm 

dEp 



(8) 



Here T(e) and /(e) = l/[e(""'^'''/'=«^ + l] are evaluated at 
thermal equilibrium. If the system is in an open circuit, 
the electrochemical potential (AV^) will be established in 
response to a temperature gradient; this electrochemical 
potential is known as the Seebeck voltage. The Seebeck 
coefficient is defined as S' = AF/AT = --C12/ Cn, where 
£11 denotes the electrical conductance Gg. 



III. RESULTS AND DISCUSSIONS 

A. Current rectification and NDC 

To calculate the tunneling current of SCQD, we adopt 
the following physical parameters: Ui = SOFq, U12 = 
IOFq, and Ei = Ep — lOFg, where Fq is a convenient 
energy unit. To simplify analysis, we ignore the magni- 
tude fluctuation of on-site Coulomb interactions and take 



Ut^UQ = 30Fo. The energy level of dot 2 {E2) is tun- 
able. 77i(2)eAya is employed to describe the energy shift 
arising from the applied voltage AVa across the junction.^ 
That means that Eg is replaced by Eg + rjfeAVa, assum- 
ing the right lead is grounded. Although the factor rji 
depends on the QD shape, material dielectric constant, 
and location, we assume that r]g is determined by the 
QD location, that is, r]i = Li/L, where Li is distance 
between the grounded electrode and the ^th QD, and L 
is the separation between the left electrode and the right 
electrode. We assume 771 = 0.6 and 772 = 0.4. This en- 
ergy level shift arising from the applied bias as observed 
in the experiment of Ref. 2 has been ignored by most 
theoretical studies of SCQD.^-i'^ 

When the value of orbital offset AE = E1—E2 is taken 
to be U2 — U12 at eAVa = 0, this would satisfy a resonant 
tunneling condition through the spin-singlet channel [the 
second channel of Eq. (3)]. It is very difficult to set up 
SCQD in the resonant condition of Ei + U12 = E2 + U2 
from experimental point of view,^ because fluctuation of 
QD size and uncertainty of its location arc hard to avoid 
in the self-assembled semiconductor SCQD.^ To reveal 
the behavior of SCQD in the off-resonance condition, we 
show the tunneling current as a function of applied bias 
for different strengths of interdot Coulomb interactions 
in Fig. 1. The dashed lines with triangle marks are 
calculated by including only the resonant function for 
the resonant channel, i.e. 



A12 
A21 



P2 



12 



ti 



|(/il-C/l2)(/i2-C/2)-i?2P' 
P5 



''|(/i2-C/2)(Ml-f/l2)-i?2P^ 

This resonant function has two poles 



(9) 



P £1 + £2 + »r ± V(e2 - ei + ^AF)^ + Uj^ 



where €1 ^ Ei + U12 + TyieAK, ^2 = E2 + U2 + ?72eAK, 
F = (Fl + Ffl)/2, and AF = (F^ - Tr)/2. Under 
the resonant condition E2 + U2 = Ei + U12, we have 
E± = Ep+ir±ti2 at zero bias {AVa = 0). Bonding and 
antibonding states are formed due to the coupling of QDs 
(ti2). We see that the dashed line matches very well with 
the black solid line (obtained with the full calculation) in 
small bias regime. This indicates that the tunneling cur- 
rent is mainly contributed by the spin-singlet resonant 
channel while the spin-triplet channel is fully suppressed. 
However, we find appreciable leakage current (J/) at high 
bias due to contribution through other channels. For 
U12 7^ IOFq, the tunneling current in the reverse bias 
is seriously suppressed. The maximum current in the 
forward bias is shifted to higher bias. Such behaviors 
are attributed to the fact that e = Ei + U12 is below 
e ~ E2 + U2 in the absence of applied bias and for neg- 
ligible ti2. The agreement between the dashed line with 
triangle marks and the solid line becomes worse when 
U12 = 6F0. The results of Fig. 1 imply that deviating 
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from the resonant condition of Ei + U12 = E2 + U2 will 
suppress the Pauli spin blockade for electron transport in 
SCQD. which plays a significant role in the application 
of spin filter. ^^^^^ When the resonant condition is met, 
the maximum current at a voltage marked by Vn^max (at 
reverse bias) is larger than that marked by Vp^max (at 
forward bias), their ratio is JR^max/ JF,max ~ 2. This is 
in quantitative agreement with experimental results re- 
ported in Ref. 2. Note that the ratio between such two 
maximum currents was reported and analyzed in Ref. 1 1 . 

Figure 2 shows the tunneling current as a function of 
applied bias for different temperatures with (a) Ei = 
E2 = Ef) = Ep — IOFq (symmetric case) and (b) 
El — E2 = U2 — U12 (asymmetric case). From Fig. 
2, we see that the current rectification of SCQD arises 
from the symmetry breaking of the carrier transport 
process due to spin blockade. This mechanism was al- 
ready illustrated using the master equation method in 
Ref. 11. We found that the ratio of reverse to for- 
ward bias maximum current {JR^max/ Jp.max) depends 
on temperature. JR^max/ Jp^max is 2.1, 1.89, and 1.74 for 
ksT = IF, ksT = 1.5Fo and ksT = 2Fo, respectively, 
and the maximum forward (reverse) current occurs at 
eAK = 2.8ro, 3.4Fo, andS.GFo (eAK = -4.4Fo, -SFq, 
and — 5.4Fo). So far, such temperature-dependent cur- 
rent rectification effects have not been reported experi- 
mentally and theoretically. 

In Fig. 2, we also notice a negative differential con- 
ductance (NDC) behavior. NDC occurs when the ap- 
plied bias is larger than Vp.max {VR.max) m the forward 
(reverse) bias regime. This NDC is attributed to the off- 
resonance behavior of QD energy levels, which can be 
tuned by the applied bias (eAV^). Therefore, it is ex- 
pected that NDC can still be observed in the absence 
of interdot Coulomb interactions. Such a NDC behav- 
ior is similar to the case of serially coupled quantum 
well, but different from the case with interdot Coulomb 
interactions. ^'''^^ The coherent tunneling current is al- 
most insensitive to temperature in high bias regime, 
which leads to a nonthermal broadening effect of the tun- 
neling current. In general, inelastic assisting tunneling 
due to phonons should also be considered for a full anal- 
ysis of the temperature or bias dependence. "'^'^It is worth 
noting that the NDC behavior of SCQD was also theoreti- 
cally studied by several workers. The NDC behavior 
of SCQD shown in Fig. 2 can be explained by the align- 
ment of the dot energy levels'^^: the current is high when 
the energy levels of different dots are aligned, but is low 
when the alignment is off. In Ref. 33, the authors also 
clearly illustrated how their NDC mechanism is different 
from those proposed in Refs. [29-32]. 



B. Charge ratchet effect 

So far, we have studied the charge transport properties 
of a single SCQD. However, it is necessary to consider the 
multiple SCQDs to achieve high spin current^"* or to cre- 



ate spin entanglement current^ in the spin filter appli- 
cation. Therefore, the proximity effect between SCQDs 
arising from the inter-dot hopping and electron Coulomb 
interactions should be included. To derive the resonant 
function (Aej) of multiple SCQDs in general based on 
the equation-of-motion method would be quite compli- 
cated. However, in the weak interdot hopping limit, we 
can apply our previous work^^ to construct the resonant 
function, Atj of multiple SCQDs by considering interdot 
Coulomb interactions for all dots,^^ while keeping the in- 
terdot hopping only between levels £ and j . When a third 
dot j' (or a charge trap impurity state) is included, the 
resonant function for the 3-dot system can be written as 

■^e.A^) = tlji^j' + bf + Cj')Y,Pm/\^m\^; £ + 3+ f- 

m 

Here the operators aji,bji and Cj> are defined differ- 
ently as in our previous work,^^ since we are consid- 
ering the effects of adding a level rather than remov- 
ing a level as in Ref. 28. Operator aj> acting on the 
terms that follow would introduce a multiplication factor 
fljv = 1 — {rij' a-) — {rij/^ff) + Cji and leave the denomi- 
nator unchanged (corresponding to adding an empty dot 
f). Operator bji would introduce a multiplication factor 
bj' ~ bj'cr + bj'ff and replace /i^ and Hj in the denominator 
by fi£ + Uij' and fJ-j + Ujji, respectively (corresponding to 
adding a singly occupied dot j'). Operator Cji would in- 
troduce a multiplication factor Cji and replace fj.e and 
in the denominator by fj.e + 2Uej' and + 2Uj.j', respec- 
tively (corresponding to adding a doubly occupied dot 
j'). Similarly, the effect of adding another dot j can be 
obtained by introducing another operator {ajn+bju+Cjn) 
to obtain the expression of j (e) for the 4-dot system. 
The procedure can be repeated for adding arbitrary num- 
ber of dots. 

Figure 3 shows the tunneling current of two parallel 
SCQDs as a function of applied bias at fesT = IFq for 
energy levels with Ei = Eq = Ep - IOFq; £ = 1,2,3,4 
(symmetric case), rji = 773 and ?/2 = '74 are 0.6 and 0.4, 
respectively. The first SCQD consists of dots 1 and 2. 
The second SCQD consists of dots 3 and 4 with dot 3 
adjacent to dot 2, while dot 4 adjacent to dot 1. There- 
fore, the hopping terms iis and ^24 are ignored. The 
intcr-SCQD electron Coulomb interactions are turned off 
in Fig. 3. Because space symmetry is maintained, the 
current spectrum is symmetrical. Comparing with Fig. 
2(a), we found that there is an extra peak in the for- 
ward (reverse) bias labeled by Jp2.max{JR2.max) result- 
ing from the electron tunneling between SCQDs (labeled 
Ji4 and ^23)- The various contributions to the total tun- 
neling current of Fig. 3(a) are shown in Figs. 3(b) and 
3(c). Fig. 3(b) shows that there are two resonant chan- 
nels with El + U12 = £4 + Uzi and ei = £4 + 2J734, 
where £1 ~ Ei + rjicAVa- The second peak at the 
bias(eAT4 = 50Fo) results from electron in state Ei tun- 
neling to state E4 + 2C/34. From the electron occupation 
numbers {Ni = N3 and N2 = N4) shown in Fig. 3(d), 
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we see that the probabihties of two electrons in E3 and 
empty E4 are high under forward bias; therefore, the 
resonant channel of £4 + 2C/34 is yielded. Because the 
SCQD is symmetric, the resonant channels under reverse 
bias are the same as those under forward bias. 

Figure 4 shows the tunneling current of two parallel 
SCQDs as a function of applied bias at fcsT = ITq for 
the asymmetric case {Ei — E2 = U2 — Ui2). Unlike the 
symmetrical case of Fig. 3, the SCQD is in the spin 
singlet state. The maximum current labeled by JFi,max 
and Jm^max are not only from intra-SCQD, but also from 
intcr-SCQD channels. These two separate contributions 
arc illustrated in Figs. 4(b) and 4(c). The resonant chan- 
nels at the conditions ei + Ui = £4 + C/4 + 2U34 and 
£4 + C/4 = £1 + 2C/i2 correspond to the current maxima 
labeled by JF2,max and JB.2,max, respectively. We note 
that the ratio of JR2,max/ JF2,max (due to inter-SCQD 
tunneling) is much larger than the ratio Jm^max/ JFi,max 
(due to intra-SCQD tunneling). This is attributed to the 
direction-dependent probability factors, which are deter- 
mined by the occupation numbers shown in Fig. 4(d). 
Highly asymmetric behavior of these occupation num- 
bers is noticed. A^i ~ is almost empty under high re- 
versed bias, which leads to a large probability weight for 
electrons entering level E4 {E2) and tunneling through 
level El {E3). So far, we have not taken into account 
the intcr-SCQD Coulomb interactions in Figs. (3) and 
(4). In a realistic system with two parallel SCQDs, the 
intcr-SCQD electron Coulomb interactions will signifi- 
cantly influence the current spectrum. Figure 5 shows 
the effects due to the presence of intcr-SCQD Coulomb 
interactions. Here, we consider U13 = [/24 = SFg and 
C/14 = U23 = 3ro. The current rectification in the 
low bias regime (cAVq < 12ro) which exists for iso- 
lated SCQD is now completely washed out, because the 
resonant condition of Ei + U12 = E2 + U2 no longer 
holds. In addition, the current spectrum is seriously 
suppressed under the reverse bias. Now, the maximum 
currents JFi,max and JF2,max result from the channels 
ei + Ui2 + Uu = £i + Ui + U24: and £i + Ui = Ei + Ui + U^i. 
The occupation numbers are shown in Fig. 5(d), which 
are useful for the analysis of charge state in each QD. The 
results of Fig. 5 imply that to control the spin charge con- 
figuration of multiple SCQDs, the proximity effect should 
be carefully taken into account. 

Recently, Si SCQDs attract serious attention for quan- 
tum bit applications due to their small nuclear-electron 
spin interaction. ^^'^^ Although Si SCQDs may have a 
longer spin relaxation time, two issues need to be ad- 
dressed: (i) The multi-valleyed nature of the Si conduc- 
tion band leads to several closely spaced energy levels 
in a Si QD. For a spherical Si QD, the six degener- 
ate valleys can be mixed by the confining potential to 
form Ai-symmetry (1 fold), T2-symmetry (3-fold) and 
E-symmetry (2-fold) states. The degeneracies may be 
further lifted by any deviation from the spherical shape, 
(ii) The presence of defect charge trap states in the ox- 
ide surrounding of Si QDs. In order to understand how 



these two issues influence the current rectification, one 
needs to calculate the resonant function of SCQD with 
closely spaced multiple energy levels. Here, we consider a 
SCQD with two energy levels per dot [Ei and E^ in dot 
1 and E2 and i?4 in dot 2). Figure 6 shows the tunneling 
current at fc^T ~ IFq as a function of applied bias for 
various values of of £"3 . When E^ and E^ are above and 
far away from EF^ the current rectification of SCQD is 
not affected [as shown in 6(a)], because these levels are 
unoccupied for all applied voltages considered. When 
i?3 is within a couple of Fg from Ep^ there exist sev- 
eral peaks in the high bias regime as labeled by JF2,max, 
jR2,max, and JF3,max- Notc that E^ IS Still above Ef, 
and the ratio JRi^max/ JFi,max is changed only slightly. 
Meanwhile, the voltages corresponding to the maximum 
current. JFi,max and Jiii,max are also nearly unchanged 
and stay close to 2ro and — 2.5ro. When E3 = Ep — Fq, 
a shoulder labeled JR2,max appears near the first peak 
jRi,max under reverse bias. The voltage corresponding 
to maximum current, JRi^max is shifted to high voltage 
(— 3Fo) as a result of the charge ratchet effect due to lev- 
els 3 and 4 (the asymmetrical behavior of A^3 and A^4 
not shown here). The resonant channels of the current 
spectrum shown in Fig. 6(c) can be analyzed as follows: 
The resonant peaks labeled JFi,max and JRi.max arise 
from electrons tunneling between level Ei and i?2. Elec- 
trons tunneling between levels E3 and E4 give rise to the 
peak labeled by JF3,max- Peaks labeled by JF2,max and 
jR2,max arise from electron tunneling between levels E3 
and E2 ■ The analysis of Fig. 6 infers that the Pauli spin 
blockade condition of the Si SCQD^^'^^ and the spin en- 
tanglement current of triple quantum dots^^ may not be 
readily implemented due to the fiuctuations of electron 
Coulomb interactions and energy levels in each QD. 



C. Thermoelectric properties 

In this section, we examine the effect of spin block- 
ade on the thermoelectric properties of SCQDs. To re- 
alize the resonant condition in the Pauli spin blockade 
regime in a SCQD junction can be challenging due to size 
fluctuation of QDs and uncertain distances between the 
fabricated QDs. In general, gate electrodes are used to 
tune the energy level of each QD to help realize the reso- 
nant condition. Figure 7 shows (a) the electrical conduc- 
tance (Ge) and (b) Secbeck coefficient (5*) of SCQD with 
Ei= Ef- lOFo and E2=Ef + 20Fo - eVg as functions 
of gate voltage and temperature. The gate voltage is ap- 
plied only to dot B. In Fig. 7(a), the four gate voltages 
Vgi, yg2 Vg3, and Vg4 tunc the energy level of dot B {E2) 
to EF^ Ef — Ui2, EF — U2 and EF — U2 — U12, respectively 
to match different resonant channels for the electron en- 
tering through El + U12. The temperature dependence 
of the peaks at Vg2 and Vgs displays a nonthermal broad- 
ening effect on the electrical conductance. These two 
peaks correspond to the resonances for the spin-triplet 
(at eVg = 30To) spin-singlet (at eVg = SOFq) states, re- 
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spcctively. It is worth noting that the magnitude of Ge 
is much smaher than 2e^//i, which manifests the effect of 
electron Coulomb interactions. In the absence of electron 
Coulomb interactions we would have Ge = 2e^//i, be- 
cause the sum of probability weights of resonant channels 
would equal to one. The peak heights of these resonances 
decrease with increasing temperature, while the widths 
are not sensitive to temperature. This behavior was first 
reported in the tunneling current measurement of SCQD 
in Refs. 5 and 6. Theoretical work for the nonthermal 
broadening of Ge was previously investigated only for 
the noninteracting system. '^^ A simple explanation for 
the nonthermal broadening effect is that the broadening 
of the tunneling current under the double-resonance con- 
dition is determined by the smaller-scale tunneling rates 
ri_2 and intcrdot coupling ti2, therefore is not sensitive to 
the larger-scale temperature variation in the distribution 
functions fL,R- The nonthermal broadening behavior of 
Ge can be useful in the application for low-tcmpcraturc 
current filtering. Our results demonstrate that a SCQD 
can function as a spin filter and low-temperature current 
filter at the same time. The two tiny peaks labeled by Vgi 
and Vgi may not be resolved in the measurement of elec- 
trical conductance Ge. However, they can be resolved in 
the measurement of Seebeck coefficient 5" as illustrated 
in Fig. 7(b). In addition, the Seebeck coefficient also 
shows a sign change with respect to applied gate voltage 
near the resonances, which arises from the bipolar effect, 
as shall be explained in details below. The inhomoge- 
neous shape of S is very different from the symmetrical 
sawtooth shape for a single metallic QD.^^ 

In Ref. 24, the electrical conductance Ge and Seebeck 
coefficient 5 of a single QD with two energy levels were 
theoretically investigated. In the case of two identical 
quantum dots coupled together, Eq. (2) can also be used 
to reveal the behavior of Ge and 5*. Figure 8 shows (a) 
Ge and (b) 5 of a SCQD with Ei = E2 = Ep - IOFq 
(symmetric case) as functions of gate voltage (applied 
to both dots) and temperature. In Fig. 8(a), the four 
gate voltages Vgi, Vg2 V^s, and Vg4 tune the energy lev- 
els of both dots {El = E2) to Ef, Ef - C/12, Ep - U2 
and Ep — U2 — f/12, while the SCQD is filled with one, 
two (quantum dot "helium" case), three, and four elec- 
trons, respectively. These peaks become broadened with 
increasing temperature, unlike in the asymmetric case of 
Fig. 7. The results of Fig. 8(a) are similar to the typical 
thermal broadening behavior of a single dot with multi- 
ple energy levels. It is noticed that S goes through zero 
when Ge reaches a maximum or minimum (which occurs 
midway between two Ge peaks). The positions of zero 
S are not affected by the temperature variation. Such a 
behavior is different from that of Fig. 7(b). In addition, 
the shape of S is also different from the sawtooth shape 
of metallic QDs with homogenous electron Coulomb in- 
teractions. We note that Ge and S in the symmetric case 
(Fig. 8) are larger than that of asymmetric case (Fig. 7). 
The above analysis should be useful for the optimization 
of the figure of merit for thermoelectric property. ^'^'^^ 



Figure 9 shows the electrical conductance Ge and 
Seebeck coefficient 5* as a function of temperature for 
two arrangements of SCQDs: (a) asymmetric case with 
El — E2 = U2 — U12 (sohd lines) and (b) symmetric case 
with El = E2 = Ep — lOFo (dashed lines). Based on Eq. 
(3), these two cases denote the spin singlet [inset of Fig. 
9(a)] and triplet states [inset of Fig. 9(b)], respectively. 
It is not easy to distinguish the singlet and triplet states 
from the electrical conductance as a function of temper- 
ature, because the electrical conductance difference in 
singlet and triplet states is small. However, the Seebeck 
coefficient provides an easy means to distinguish the spin 
singlet and triplet states, since the sign of Seebeck coef- 
ficient changes for the spin triplet state as temperature 
increases, but for the spin singlet state, the Seebeck co- 
efficient is always negative, because the electrons from 
the left electrode (hot side) diffuse into the right elec- 
trode (cold side) through the resonant channels above 
the Fermi energy, which leads to the built-up of negative 
Ay in order to reach the J = condition for open circuit. 
Note that the resonant channel of Ei + U12 = E2 + U2 
has no contribution to the Seebeck coefficient (£12 — 0). 
For the spin triplet state, zero Seebeck coefficient occurs 
at To, i.e. S{Tq) = 0, which indicates that the current 
arising from temperature gradient can be self-balanced 
without electrochemical potential (Ay = 0). On the 
other hand, the Seebeck coefficient is positive when the 
holes from the left electrode (hot side) diffuse into the 
right electrode (cold side) via the resonant channels be- 
low Ep- Here, we have defined the unoccupied states 
below Ep as holes. Consequently, the sign change in S 
is attributed to the competition between tunneling cur- 
rents due to electrons and holes (so-called bipolar effect). 
This implies that we can control the current direction by 
manipulating the equilibrium temperature. In addition, 
we can distinguish the spin singlet and triplet states by 
measuring the temperature dependence of the Seebeck 
coefficient. 



IV. SUMMARY 

We have used a two-level Anderson model to describe 
a SCQD connected to metallic electrodes and calculated 
the tunneling current within the framework of nonequi- 
librium Green's function technique. In the Coulomb 
blockade regime, we have derived a closed-form expres- 
sion for tunneling current. In the nonlinear response 
regime, the electron spin singlet and triplet states can be 
distinguished by the current rectification behavior aris- 
ing from the space symmetry breaking, whereas such a 
current rectification effect in the spin blockade process 
is suppressed with increasing temperature. We have also 
studied the proximity effect between two parallel SCQDs 
and analyzed how the charge trapping states influence 
the current rectification of SCQD. In the linear response 
regime, the electrical conductance Ge and Seebeck co- 
efficient S are analyzed. It is not easy to distinguish 
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Ge due to transport through the spin singlet or triplet 
states from its temperature behavior. However, the tem- 
perature dependence of Seebeck coefficient can clearly 
reveal the spin configuration by examining the sign of 
S. For example, we observe a sign change in S for elec- 
tron transport through the spin triplet states, while S 
is always negative for transport through the spin singlet 
state. Thus, we conclude that the measurement of See- 
beck coefficient provides a much better means for resolv- 
ing the resonant channels than measuring the electrical 
conductance. In addition, we see a nonthermal broad- 
ening effect of tunneling current for the spin blockade 
process. This indicates that SCQD can simultaneously 
act as the spin filter and low-temperature current filter. 
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Appendix A. 

The expression of the tunneling current through seri- 
ally coupled quantum dots (SCQD) is derived based on 
the seminal work by Meir and Wingreen ' . and Jf{ 
denote, respectively, the tunneling current of electrons 
leaving the left and right electrodes, which are expressed 
by 



Jl = 
and 



J, 



B. 



J de^,{e)[2h{e)lmGl^,ie)-^G<,{e)], (A.l) 



der.2{e)[2Me)lmGl^{e) - tG<,{e)]. (A.2) 



According to expressions of Eqs. (A.l) and (A.2), tun- 
neling current is determined by the on-site retarded 
and lesser Green's functions (GFs). Note that we con- 
sider the SCQD system to be spin degenerate, thus 
have suppressed the spin index for convenience. Here 
G^ j{t) = —id{t){{d(,a-{t), d'j ^}) and its Fourier transform 

one-particle GF. 
It is nontrivial to obtain an exact solution of G^ ^ (e) for 
the system Hamiltonian when the coupling with leads is 
present.^* We solve the on-site retarded and lesser GFs in 
the Coulomb blockade regime. One-particle GF G^ j{e) 
can be obtained by solving a closed set of equations of 
motion. For coupled QDs, it is convenient to define vec- 
tor GFs as 

- ( StitI ); - 

(ni^niird2a\dl^) 
{n2aniadia\dl^) 

^ \ {n2,Tniad2a\dl^) 

{n2cin2aniadla\d\^) 

{n2cTnian-iad2a\d\^) 
denote the number of particles involved, and we have de- 



{nigdia\d\^) 
{nigd2a\d\„) y ' 



M2) ^ / {n2adla\d\^) 



1 



i=l,2,12, or 21, G, 



(3) 



and 



where the superscripts 



fined n(jg = ^i-dj^; i,j = 1,2. a denotes the opposite 
of spin a. 

From Eq. (1) we obtain the equations of motion that 
relate g^^) to g'-^'^'s, then to G^^^'s, and finally to g^^\ 
which self-terminates. The one-particle GFs satisfy 

HoG^^'' = ( J ) + + Wef ^ + U^2GP , (A.3) 



where we have defined Hq 



-tl2 M2 

fie^e- Ee^, + ire/2; ^ = 1, 2 (F, 27r \Vk,e\^S{e - 
fk) arising from the QD coupled to electrode), U — 

jj ^ ^ ^^'^ ^ ^ (^^0^ U ^ ' "^^^ two-particle 
GFs satisfy 

{Ho-U)G['' = ^(J" ) +h2{G[l'-Gif)+UG^'' + U^2G['\ 

(A.4) 

(A.5) 

{Ho - u^2)Gr = ( ) +ug? +uGf\ (a.6) 



,(2) / N2a 



{no-U)G2 - I 

,{2) / N2cr 
,(2) _ / (?ll2a) 



tl2iG[V -GiV)+UG^'^ +Ul2Gi'\ 



{Ho-U-Ae)Gr2 = 







+ ti2{G[^^ ~Gi''^)+hl'G[l\ 
(A.7) 

(% -w + Ae)g(f = ( ^""'^^^ ] -ti2(er^ -af ) 



where Ae = E2 ~ Ei and W = 
The three-particle GFs satisfy 

(Ho - W-t/l2)ef^ -- ' 



{Ho~U~Ui2)G 



(3) 



-Ci2(t,1ct 
C2 

'Cl2cr,2a 



(A.8) 

72 \ 

-UiJ- 

+ UG^*\ (A.9) 
+ VIG^*\ (A.IO) 



(% - W - W)af> = ( ""^'^'^ ] + Ui2G^^\ (A.ll) 



fj,i - Ae -U2 - U12 -ti2 \ ^(3) 

-ti2 fl2 - Ae - U2 + Ul - 2f/i2 



G12 



C2.,,2. ) +ti2(ef^-Gf ^ 



Ml + Ae - {/i -f t/2 - 2(7i2 



-ii2 



-t-L2 

fi2 + Ae-Ui- Ur. 



C2tT,21& 
-Ci2ct,21ct 



-ti2(Gr' -ef ), 



(A.12) 



,{3) 



(A.13) 



where we have defined the c-numbers ct = {ni_„ni^g), 
and Cia,i'a' = {niarii'a')] hi' = 1,2,12,21. Finally, the 
four-particle GFs satisfy 

^ ^ V -\n12an2ania) J ^ ' 
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('21 

From Eqs. (A.4) and (A.5), we see that the terms 
and Q^'i only give a small correction of order t\2/U to 

(2) (2) 

Ql and Q2 , while those off-diagonal expectation values 
such as {ni2a) that are of the first order in ii2 in Eq. 
(A. 6) give a correction of tyilV; thus, these terms can be 
ignored in the limit of small tyijU . If we further make 
the approximation {n2an2gnig) ~ C2Nig (valid again in 
weak interdot coupling case), then the one-particle GFs 
can be written in simple closed forms. 

Ftom Eqs. (A. 3-14), we found that the one-particle re- 
tarded GF can be simply decomposed into a sum of con- 
tributions from eight spin-charge configurations of the 
SCQD system in the condition of t\2llJi ^ 1 (after ig- 
noring those small contributions proportional to ti2/U .). 
It reads 



7(1) 



= -H^^Vi + (% - U)-^V2 + {-Ho - Ui2)-^Vz 
+ (Ho -U- UriY^Vi + (7^0 - Uy^Vr, 
+ {Ho -U- Uy^Pe + {Ho-U- C/i2)~^P7 
+ {Ho-li -U -U^y^Vs, (A.15) 



where Vn 



Pm 





; Pin's are probability weights de- 



fined in Eq. (3). The intradot one-particle and two- 
particle retarded GFs are given by 



Pi 



P2 



and 



I^Le - tl^lixj {fii - Ui^j) - tl^l (mj - Uj) 

P3 

- Utj)-tl2l{iJ.j - Uj^l) 

Pi 

Pb 

P6 

(/if -Ui- Uij) - tl^/ipLj - Uj - Uj^e) 

P7 

P8 

(A.16) 



P5 



Pa 



(M£ -Ue- Uej) - tl^/ijij - Uj - Ujj) 

+ El 

im - Ue - Ui,y - tl^l {^Lj - 2Uj.i) 

+ 

ifie -Ut- 2Ut,j) ~ tl^/if^j - Uj - 2Uj,e) ' 

(A.17) 

The index j denotes the jth QD, and j 7^ £ in Eqs. (A.16) 
and (A.17). 



Here we shall derive the one-particle lesser GFs Gfi 
and G22 required for the current formulae (jA.llA.2p 
directly from the nonequilibrium equation-of-motion 

method. Define Qfj{t) = i{d}j ,^{Qi)di^a{t)) and its Fourier 
transform as {di^a\d?^ g.)^; similarly for higher lesser GFs. 

Using the nonequilibrium equation-of-motion method"^^ 
and treating the coupling to the electrodes to lowest or- 
jgj,8,28^ we can readily obtain a closed set of equations of 
motion: The one-particle lesser GFs satisfy 

Hoe<<'' = E<e;''(^' + wg<^'' +ug<'^^'^ + UuQt''^\ (a.is) 



where we define in 2 x 2 matrix form 







E< 



with = iTa{e)fa{() for a = 1,2 and one-particle 
advanced GF vectors in Eq. (A.18). The two- 

particle lesser GFs satisfy 



.<^a(2) . , ,^<(3) . rr „^<{3) 



.<ga(2) 

The three-particle lesser GFs satisfy 



{-Ho - Ui2)gs'-^'' = E^g;'^' +ugy" +ug^ 



+ u^g; 

<(3) I 77^<{3) 



(A.19) 

(A.20) 
(A.21) 



(•Ho 


-U 




= E<gf'^ +ug^'-*\ 


(A.22) 


(•Ho 
(•Ho 


-U 
-U 


- u,2)g^^'~> 
-u)gf'^^ 


= ^<gf'^ +ug^^^\ 
E<g3"(»' + [/i2g<w. 


(A.23) 
(A.24) 



The four-particle lesser GFs satisfy 



(A.25) 



The approximations made here conform to those at de- 
riving the retarded GFs in the weak ti2 limit. Similarly, 
we found that the lesser GFs can equally be decomposed 
into a sum of contributions. It reads 

Ho^Y.<H*o-^Vi + {no-uy^j:<{n*o-uy^P2 
+{'Ho - Ui2y^j:<{H*o - Ui2y^r3 
+{no-u-Ui2y'^<{n*„-u-Ui2y'Vi 
+{no-uy'^<{n*o-uy'r5 

+{nQ-u -uy^j:<{n*o-u -uy^Pe 
+inn-u -Ui2y'^<{K-u -Ui2y^V7 

+{no -u-u- Ui2y^^^{n*o -u-u- Ux2y^Vs. 

(A.26) 

Straightforward algebra leads to 

= !]< 



Pi 



P2 



\^Ji, - i?2/M,f - u,,,) - tl2/{ii, - u,y 

P3 



\i^le-Uf.,)-tl2/{p,-U,Jy 

P4 

'\i^ie-2Ue,J)-tl2/i^^J-U,-UJ^,) 
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P5 



\{fie-Ue)-tl,/{^,,-U,.eW 

P6 

P7_ 

■ |(W -Ue- C/,,j) - tlJ{^ij - 2C/,- ,)|2 



-Ue- 2Ue^,) - tlJip., - U, - 2Uj^eW 



+^12 
+ 



Pi 



P2 



P3 



|(w-C/£j)(Mi-C/j,£)-i?2P 
Pi 

^{^ie-2Ue^,){^i,-U,-U,J)-tl^\^ 

P5 

^{^ie-Ue){^l,-U,.e)-tl^\^ 

P6 

'\{tie-Ue- Ue,j)itxj - U, - U,.e) - tf^l^ 

P7 

'\i^^e-Ue-Ue,,){p,-2U,,e)-tl^\' 
Ps 



\ifie -Ue- 2Ue,j)ifij - U, - 2U,,e) -tl^? 



(A.27) 



The index j in Eq. (A.27) denotes the jth QD, and j ^ 
Inserting the diagonal GFs defined by Eas. (|A.16IA.27|) 
into Eqs. (jA.l|A.2p yields 



Jl = ^ y de[fL{e) - /i?(e)]ri(e)r2(e)A2(e), (A.28) 
and 

jR = lj de[fR{e) - /L(e)]ri(e)r2(6)^2i(e). (A.29) 

Furthermore, using J = Jl = {Jl — Jr.)/2 and = 
— Jr, we can symmetrize the current as 



J 



2e 



der{e)[h{e) - fnie)] 



(A.30) 



where we define the transmission factor 

r(e) = ri(e)r2(e)(A2(e)+ Ai(e))/2. (A.31) 



The factor 2 in Ea. (|A.30|) accounts for the electron spin 
degree of freedom. An and A21 in Eq. (A.31) are called 
resonant functions, which arc given by 



= (/i£ - Ue - Uej){nj - 2Uj,e) - tj^, 
Us - if^e -Ue- 2Ue,j){fij - U, - 2Uj,e) - tj^- 
We emphasize that the current formula of Eq. (A.30) is 
correct up to the second order in ti2 . Taking into account 
those off-diagonal expectation values such as {ni2a) can 
give a more accurate current formula. However, in the 
limit of our interest where ti2 is much smaller than the 
other energy scales, we do not expect any major differ- 
ence. 

In this equation-of-motion framework, the expectation 
values (ne^a) and ce can be readily computed via {ne^a) = 
J deGfg/2TT and ce = J deG^g f /2tt. Now the thermally 
averaged occupation number is given by 



No 



{ne..) = / ^G<,(e) 



Ce 



^ " / f /Ke)IniG?,,(e); (ti2«r,),(A.33) 



de 



feie)lmGf',{e); {t^ « T,). (A.34) 



Finally, we compare our results with those in Ref. 17, 
where Yuan et al found the current formula in the small 
ti2 limit to be 

(A.35) 

From this current expression we understand that the den- 
sity of states for the SCQD is determined by the product 
of the density of states of each QD, as expected in the 
small ti2 limit. Their solution for the retarded GF is 



{e-Ee..){e-E,,,-Ue) 

; h li ( 2 

e-Ee- Ue{l - Ne^-a) ' 

1 - Ne,-a , Ne.-a 



e - Ee,a + ^Ff /2 e - -Ue + iVel2 ' 

for Ue 3> Ff. It can be readily checked that our result 
is identical to Yuan et aFs work if we turn off interdot 
Coulomb interactions. 



Aejie)=t 



12 X! 



(A.32) 



where the denominators for the eight configurations are: 

Hi = i^ie^j - tj^, n2 = (ue - Ue^j){iij - Uj) - tj^, 

Hs = {jie - Ue,j){fij - Uj,e) - tf.^, 

Hi = (lie - 2Ue,j){H - - Uj^e) - ^2, 

Hs = {fie - Ue){fij - U^j) - tj^, 

He = ifJ-e - Ue - Ue,j){fJ,j - Uj - Uj^e) - t^^, 



10 



^ E-mail address: nitkuo@ee.ncu.edu.tw 



* E-mail address: yiachang@gate.sinica.edu.tw 



^ H. Haug and A. P. Jauho, Quantum Kinetics in Trans- 
port and Optics of Semiconductors (Springer, Heidelberg, 
1996). 

^ K. Ono, D. G. Austing, Y. Tokura and S. Tarucha, science 

297, 1313 (2002). 
^ K. Ono and S. Tarucha, Phys. Rev. Lett. 92, 256803 

(2004). 

* A. C. Johnson, J. R. Petta, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Phys. Rev. B 72, 165308 (2005). 

5 N. G. van der Varrt, S.F. Godijn, Y. V. Nazarov, G. J. P. 

M. Harmans,J. E. Mooij, L. W. Molenkamp,G. T. Foxon, 

Phys. Rev. Lett. 74, 4702 (1995). 
® W. G. v. Wiel, S. D. Franceschi, J. M. Elzerman, T. Fuji- 

sawa, S. Tarucha, and L. Kouwenhoven, Rev. Mod. Phys. 

75, 1 (2003). 

Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

* B. R. Bulka and T. Kostyrko, Phys. Rev. B 70, 205333 
(2004). 

^ P. Pals and A. Mackinnon, J. Phys.: Gondens. Matter 8, 
5401 (1996). 

^° S. A. Gurvitz and Y. S. Prager, Phys. Rev. B 53, 15932 
(1996). 

J. Fransson, and M. Rasander, Phys. Rev. B 73, 205333 

(2006) . 

^2 B. Muralidharan and S. Datta, Phys. Rev. B 76, 035432 

(2007) . 

" J. Inarrea, G. Platero, A. H. MacDonald, Phys. Rev. B 76, 
085329 (2007). 

" Q. F. Sun, Y. X. Xing, and S. Q. Shen, Phys. Rev. B 77, 
195313 (2008). 

R. Hornberger, S. KoUer, G. Begemann, A. Donarini, and 
M. Grifoni, Phys. Rev. B 77, 245313 (2008). 
^® P. Trocha, I. Weymann, and J. Barnas, Phys. Rev. B 80, 
165333 (2009). 

" R. Y. Yuan, R. Z. Wang and H. Yan, J. Phys. Gondens, 

Matter 19, 376215 (2007). 
1^ W. Gong, Y. Zheng, Y. Liu, and T. Lu, Phys. Rev. B 73, 

245329 (2006). 

1^ D. M. T. Kuo, S. Y. Shiau, and Y. G. Ghang, 

|arXiv:1101.5f51 
^° A. Vidan, R. M. Westervelt, M. Stopa, M. Hanson, and A. 

C. Gossard, Appl, Phys. Lett. 85, 3602 (2004). 

2^ D. Culcer, L. Cywinski, Q. Z Li, X. D. Hu, and S. Das 
Sarma, Phys. Rev. B, 82, 155312 (2010). 
S. Das Sarma, X. Wang, and S. Yang, Phys. Rev. B, 83, 
235314 (2011). 

P. Murphy, S. Mukerjee, J. Morre, Phys. Rev. B 78, 161406 
(R) (2008). 

D. M. T. Kuo and Y. C. Ghang, Phys. Rev. B 81, 205321 

(2010) . 

Y. Dubi, and M. Di Ventra, Rev Modern Phys 83, 131 

(2011) . 

X. Zianni, Phys. Rev. B 78, 165327 (2008). 

R. Sanchez and M. Buttiker, Phys. Rev. B 83, 085428 

(2011). 

D. M. T. Kuo and Y. G. Chang, Phys. Rev. Lett. 99, 



086803 (2007); Y. G. Chang and D. M. T. Kuo, Phys. 
Rev. B 77, 245412 (2008). 
2^ B. Dong, H. L. Cui, and X. L. Lei, Phys. Rev. B 69, 035324 
(2004). 

B. Wunsch, M. Braun, J. Konig, and D. Pfannkuche, Phys. 
Rev. B 72, 205319 (2005). 

L. Djuric, B. Dong, and H. L. Cui, J.Appl. Phys. 99, 063710 
(2006). 

J. Aghassi, A. Thielmann, M. H. Hettler, and G. Schon, 
Phys. Rev. B 73, 195323 (2006). 

J. N. Pederson, B. Lassen, A. Wacker, and M. H. Hettler, 
Phys. Rev. B 75, 235314 (2007). 

D. S. Saraga and D. Loss, Phys. Rev. Lett. 90, 166803 
(2003). 

L. Oroszlany, A. Kormayous, J. Koltai, J. Cserti, and G. 
J. Lambert, Phys. Rev. B 76, 045318 (2007). 

C. Niu, D. L. Lin, and T. H. Lin, J. Phys. Gondens. Matter 
11, 1511(1999). 



11 



Figure Captions 

Fig. 1. Tunneling current as a function of applied bias 
for the variation of U12 at temperature ksT = IFq, and 
Tl = Tr = 0.5r. Ei^Ep- lOFo and E2 = Ep ~ SOFo. 
Jo = 2 X lO-^eVo/h). Other physical parameters = 
Tr = O.SFo, and ti2 = O.lFo. 

Fig. 2. Tunneling current as a function of applied 
bias for different temperatures, (a) Ei — E2 ~ Ep — 
lOFo (space symmetrical case), and (b) Ei ~ Ep — IOFq 
and E2 = Ep — SOFq (space symmetry breaking) . Other 
physical parameters are the same as the black line of Fig. 
1. 

Fig. 3. Tunneling current as a function of applied bias 
at ksT = IFq in the absence of inter-SCQD electron 
Coulomb interactions. ti2 = <34 = 0.025Fo and ti4 = 
h2 = O.OlFo. Tl ^ Tr ^ 0.5Fo. Ee = Ep - IOFq, 
Ue = 30Fo. and U12 = U34 IOFq. 

Fig. 4. Tunneling current as a function of applied bias 
at fcsT = IFo in the absence of inter-SCQD electron 
Coulomb interactions. Ei ^ E3 = Ep — IOFq, and E2 = 
E4 = Ef — SOTq. Other physical parameters are the same 
as Fig. 3. 

Fig. 5. Tunneling current as a function of applied bias 
at ksT = IFo in the presence of inter-SCQD electron 
Coulomb interactions. U14 = C/32 = 3Fo and U13 = 



U24 = 5Fo Other physical parameters are the same as 
Fig. 4. 

Fig. 6. Tunneling current as a function of applied bias 
for different energy levels of E3 at E4 ~ Ep + IOFq and 
ksT = IFo for the SCQD with two energy levels per dot 
{El and E3 in dot 1 and E2 and E4 in dot 2). Other 
physical parameters are the same as Fig. 5. 

Fig. 7. (a) Electrical conductance Ge and (b) Scebcck 
coefficient S of the SCQD with Ei = Ep ~ IOFq and 
E2 ~ Ep + 2OF0 — eVg as functions of gate voltage for 
various temperatures. Other physical parameters are the 
same as the black line of Fig. 1. 

Fig. 8. The electrical conductance Ge and Seebeck 
coefficient 5* as a function of gate voltage for different 
temperatures at Ei = E2 = Eq = Ep + 20Fo — eVg. 
Other physical parameters are the same as the black line 
of Fig. 1. 

Fig. 9. The electrical conductance Ge and Seebeck 
coefficient as a function of temperature for the spin 
triplet states (dashed lines) and spin singlet state (solid 
lines). For dashed lines, the SCQD is spatially symmetric 
(triplet state) with Ei = E2 = Ep — IOFq. For solid lines 
the SCQD is spatially asymmetric (singlet state) with 
Ei^ Ep- lOFo and E2^Ep- SOTq. Here U12 = IOFq. 
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